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We study the coupling of the gravitational action, which is a linear combination of the 
Hilbert-Palatini term and the quadratic torsion term, to the action of Dirac fermions. The system 
possesses local Poincare invariance and hence belongs to Poincare gauge theory with matter. The 
complete Hamiltonian analysis of the theory is carried out without gauge fixing, which leads to 
a consistent geometrical dynamics with second-class constraints and torsion. After performing 
a partial gauge fixing, all second-class constraints can be solved, and a connection dynamical 
formalism of the theory can be obtained by a canonical transformation. Hence, the techniques of 
loop quantum gravity can be employed to quantize this Poincare gauge theory with non-zero torsion. 
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I. INTRODUCTION 



O . General Relativity (GR) has been very successful in describing universe at large scales. However, it is believed that 
we have to develop a quantum theory of gravity for a consistent description of nature. One of the reasons that classical 

$H . GR cannot be consistent can be seen from the Einstein's equations which relate gravitational and matter degrees of 
freedom. While the gravitational part is classical and is encoded in the Einstein tensor, since matter interactions are 
very well described by quantum field theory, we need to use some quantum version of the stress energy tensor for the 
^-H ' matter part. This would imply that a consistent coupling of matter and gravity for all energy scales requires both of 

I> ] them to be quantized. 

. Einstein's equations can be obtained via an action principle starting from the first-order Hilbert-Palatini action. 

However, if we consider fermionic matter sources, the equations of motion from this action will not provide the torsion- 
, free condition of vacuum case. Hence, we have to either allow for torsion or make some suitable modification of the 
action. (See [l| and references therein for a comprehensive account of torsion in gravity). So, if one wants to start 
, with first-order action, it is very possible that quantum theory of gravity would incorporate torsion in its formalism in 
order to consistently couple gravity to fermions. Among various attempts to look for a quantum gravity theory, gauge 
theories of gravity are very attractive since the idea of gauge invariance has already been successful in the description 
of other fundamental interactions. Local gauge invariance is a key concept in Yang-Mills theory. Together with 
Poincare symmetry, it lays the foundation of standard model in particle physics. Localization of Poincare symmetry 
leads to Poincare Gauge Theory(PGT) of gravity. One of the key features in PGT is that, in general, gravity is not 
; ^ ■ only represented as curvature but also as torsion of space-time. GR is a special case of PGT when torsion equals zero. 
Ci , PGT provides a very convenient framework for studying theories with torsion. A number of actions which satisfy 
local Poincare symmetry have been analyzed by various researchers ( [5| provides a comprehensive review and bibliog- 
raphy of the progress made in PGT). However, one of the drawbacks of PGT is that its Hamiltonian formulation is 
usually very complicated. Although Hamiltonian analysis is performed for many models in PGT, the results are at a 
formal level without explicit expressions of the additional required second-class constraints. From the point of view 
of canonical quantization, it is essential to have a well-defined consistent Hamiltonian theory at the classical level. 
Such an ingredient is missing if we want to incorporate torsion into candidate quantum gravity models constructed 
from PGT. Moreover, the internal gauge group in PGT is in general non-compact, while most of the standard tools 
developed in quantum field theory apply to gauge theories with compact gauge groups. 

There exists a well-known SU{2) gauge theory formulation of canonical GR where the basic variables are the 
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densitized triad and Ashtekar-Barbero connection. A candidate canonical quantum gravity theory known as Loop 
Quantum Gravity (LQG) i^iS] can be constructed starting from the connection dynamical formulation. Moreover, 
LOG can also be extended to some modified gravity theories such as, f{R) theories 0, and scalar-tensor theories 
[11]. However, the action of GR from which the connection dyanimics can be derived is not the standard Hilbert- 
Palatini action. An additional term known as the Hoist term has to be added to the standard Hilbert-Palatini action 
in order to rewrite GR as a SU{2) gauge theory [l^, ES]- It is customary to multiply the additional Hoist term with a 
coupling constant 7 known as the Immirzi parameter. Classically these two actions are equivalent since the additional 
Hoist term does not affect the equations of motion although it is not a total derivative. The parameter 7 does not 
appear in the classical equations of motion. This is because the Hoist term differs from a total derivative known as 
the Nieh-Yan term 14 [ by a term quadratic in torsion (for the exact relations between them see [H, [3). Since the 
torsion term is zero when there is no fermionic matter, the Nieh-Yan term and the Hoist term are same, and hence the 
connection dynamics obtained from adding either term to the Hilbert-Palatini action would be equivalent. It has been 
shown that a SU{2) gauge theory can also be constructed from an action containing the standard Hilbert-Palatini 
term and the Nieh-Yan term [l^. However, when there are fermions, the term is not zero and the the difference in 
the Hoist term and the Nieh-Yan term shows up. In [l3l it was found that adding the standard fermion action along 
with the Hoist term leads to equations of motion which depend on 7 and are therefore not equivalent to standard 
GR with fermions. The difference arises because the Hoist term is not a total derivative. In [15| it was shown that 
there is no such issue if the full Nieh-Yan term is used. An alternative possibility of modifying the fermion action 
to be non-minimally coupled has been analyzed in detail in [isl fioj . The additional piece in fermion action cancels 
the contribution of the Hoist piece if the coupling constants are chosen accordingly. (See H^j for a recent account of 
these issues) . In the absence of direct experimental or observational evidence of quantum gravity and of torsion, it is 
not clear which action should be the appropriate starting point for quantization, particularly from the perspective of 
LQG. It is therefore very important to study all the different possibilities. However to apply the LQG techniques, it 
is essential to first reformulate these candidates as gauge theories with a compact gauge group. 

In this series of works, instead of the Hoist piece of the Nieh-yan term, we consider the piece. In [2lj we 
considered the vacuum case, i.e. an action with only this along with the standard Hilbert-Palatini term. Instead 
of the Immirzi parameter 7, we use an arbitrary coupling constant a which is related to 7 by numerical factors. There 
it was shown that, although we started from an action with explicit torsion dependence, the conditions of the closure 
of the constraint algebra imply that torsion is zero, and hence we go back to standard GR. This is consistent with the 
results that there is no torsion in the absence of spinors. However, in that case it was not possible to obtain a SU (2) 
gauge theory. The variables we choose are motivated by PGT. But unlike other analyses in PGT we obtain explicit 
expressions of the second-class constraints. 

In this paper, we add Dirac fermions to the action and apply the techniques developed in (2]| to carry out the 
Hamiltonian analysis. We consider the fermions to be non-minimally coupled, because the term is not a total 
derivative and indeed, by proper choice of the two coefficients, the contribution of the additional non-minimal piece 
is cancelled by the contribution of the torsion piece in the closure of the constraint algebra. Also the relation between 
torsion and the fermions we obtain is the same as the one obtained in [TE^ with Nieh-Yan term and minimally coupled 
fermion action. To the best of our knowledge, this is the first action with explicit torsion terms which has been 
reformulated as a Hamiltonian SU (2) gauge theory. The new connection we obtain has some properties which are 
novel and very different from the other connection dynamics formulations obtained so far. The classical system we 
obtain in this paper can subsequently be loop quantized using the tools already developed in LQG and will lead 
to a new and inequivalent quantum theory. Also, Hamiltonian formulation of theories with torsion are usually very 
complicated. We think that the techniques developed in this and the previous paper [2l| can be used for analyzing 
other similar actions with torsion terms. If that is possible, then the general programme of loop quantization can be 
applied to a much wider class of theories which include torsion. 

The paper is organized as follows. In Section pi)) we give the explicit expression of the action we start with and 
our definitions and conventions. In Section (jIII[) we perform a 3 -1- 1 decomposition of this action and perform the 
Hamiltonian analysis without fixing time gauge. Having obtained a consistent Hamiltonian system, we fix time gauge 
and then solve the second class constraints in Section (jIVp . Fixing the time gauge also breaks the 5*0(3, 1) gauge 
invariance to SU{2). Then in Section (|V| we define a new connection which is conjugate to the densitized triad to 
obtain a SU{2) gauge theory. Our analysis has several novel and peculiar features. We conclude with a discussion of 
these and some comparison of our results with those obtained using the Hoist and Nieh-Yan terms in Section (jVIl) . 
We will restrict ourselves to 4 dimensions. The Greek letters ^^v . . . refer to space-time indices while the uppercase 
Latin letters I,J... refer to the internal 5*0(3, 1) indices. Our spacetime metric signature is ( — h Later when 

we do the 3-1-1 decomposition of spacetime, we will use the lowercase Latin letters from the beginning of the alphabet 
a,b, . . . to represent the spatial indices. After we reduce the symmetry group to SU{2), the internal indices will be 
represented by lowercase Latin letters from the middle of the alphabet 
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II. THE ACTION 

In this paper we consider an action which has three pieces, a Hilbert-Palatini term, a term quadratic in torsion and 
a term for the massless fermionic matter. It reads 



S — Shp + cxSt + Sm 



(1) 



where 



Shp = J d^x ei?, = J d'xee'ie-'jRj''{LoJ'') , 

Sm = *e / A\ [A(l + ipf)-f''D^\ - d;Jj^{1 + il3-f^)X] 



Here Cj is the tetrad, e denotes the absolute value of the determinant of the co-tetrad, uj"^ is the spacetime spin- 
connection which is not torsion-free, e'^'^'"^ denotes the 4-dimensional Levi-Civita tensor density, and the covariant 
derivatives in the fermion action read. 



1 



D^X 



dnX — — X UJ 



IJ 



(TIJ. 



Note that we denote 7^ = 'j^Cj with 4-dimensional Dirac matrices 7^ , ct/j :— | [77,77] and 75 := «7o7i7273- Our 
conventions regarding the Dirac matrices and their properties are given in Appendix ([XJ. Note also that A and 
A := A^7'', representing the fermionic degrees of freedom, are 4-dimensional row and column vector respectively. 
Further, 



R 
T 



= di 



I J 

I — T ^ 



-'^[/|j|e;^] 



J 

\K ; 



(2) 
(3) 



are the definitions for curvature and torsion respectively. ^ The action ([T]) is invariant under local Poincare transfor- 
mations. We will be working in the first-order formalism and hence both the co-tetrad and the spin connection 
w^^'' are treated as independent fields. Our covariant derivative Z?^ acts in the following way: 



- UJ 



J^v ■ 



Note that the coupling parameter a in action ([1} is a non-zero real number. The parameter /3 in the matter action 
denotes nonminimal coupling and with /3 = we get back minimally coupled Fermion action. The parameter /3, in 
general, has no relation with the parameter a. However in the subsequent Hamiltonian analysis, we have to choose 
/3 = ^ for the closure of the constraint algebra. We therefore adopt that relation between the two parameters from 
here onwards. 

By variation of the action we can find the equations of motion for the various fields, namely e^, ujfj'\ X and A. In 
particular the equation of motion for the connection uo^'-' is given by 



-^{VJKVIL - VIKVJl) - <^IJKL 



iee^A757LA 



-^{VIKVJL - VJKVIl) + <^IJKL 



= 0. 



(4) 



It is easy to see that, as expected in the presence of spinors, the connection is not torsion free, i.e., the equation of 
motion do not imply D^pC^^ =0. In 21 1, the Hamiltonian analysis of the action ([1]) without the matter part was 
carried out. In that case, the Lagrangian equations of motion showed that torsion was zero on-shell although the 
action has explicit torsion terms. In the next section we will carry out a similar analysis with action ([!]) where the 
torsion is expected to be non-zero. 



^ Our conventions of symmetrization and antisymmctrization are A^""^^ := A'^^ + A*"" and A^°-^^ := A"*" — A*"* respectively 
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III. HAMILTONIAN ANALYSIS 

We shall perform the Hamiltonian analysis similar to what was done in [2lj for the action without the matter 
term. Recall that in the Hamiltonian formulation of Hilbert-Palatini theory the basic variables are the 50(3, 1) spin 
connection u)^"^ and its conjugate momentum. It is well known that this formulation contains second-class constraints. 
Since our action ([ij contains the other term which explicitly depends on torsion, we expect that there will be another 
pair of conjugate variables and the second-class constraints will be somehow different from the Hilbert-Palatini case. 
It is well known that in the absence of fermionic matter, torsion is zero. In the analysis of ^21i] . this was obtained as 
a condition for the closure of the constraint algebra. Owing to the presence of the fermion term in the action, here 
torsion will not be zero. However the conditions for closure of constraint algebra will determine how the torsion and 
the spinorial degrees of freedom are related. 



A. 3+1 Decomposition 



To seek a complete Hamiltonian analysis, we perform the 3-1-1 decomposition of our fields without breaking the 
internal 50(3, 1) symmetry and also without fixing any gauge. To identify our configuration and momentum variables 
for performing Hamiltonian analysis, we can rewrite the three pieces in the action as: 



S 



HP 



d^x 



1 



IJ 



(5) 



a d X 



Sm — d X ie 



(a(1 + i^i'WdtX - a*A7*(l + *f 7')a) + \ (a(1 + z|7')7 



+ ( A(l + z|7')7'^^aA - DaXi'^il + *|7')A 



(6) 



a 



(7) 



We can read off the momenta with respect to ujJ'\ e^, A and A respectively as 



J := ae Ubed 



.a 



H := ieX(\ ^ i^-f')e\-i^ 



H 



-^ee*7^(l-H^-7')A, 



(8) 
(9) 



where e"''^ denotes the 3-dimensional Levi-Civita tensor density, and we have used the relation 7^ = e'^. For our 
analysis we shall use a standard parametrization of the tetrad and the co-tetrad fields as in |22l|. This is the same 
parametrization used in the Hamiltonian analysis of the first two terms of our action in [2l|. The details of the 
parametrization and some useful identities are given in Appendix (|B]) . 

After some manipulation and neglecting the total derivatives, the pieces of the action can be written in this 
parametrization respectively as 



Sm — 



d^x 



d^x 



d^. 



2g ^^IK^^JL^ ^ab ^2 IJ^ab ^a^^IJ 



na^A + mix - (^nj,^ (na^'^z?,A - a^a^t^h) 

-HiV" (Hi^.A + A^H) + ic^/^ (Aa/jH - Ha/jA) j 



(10) 
(11) 



(12) 



Collecting them together we can rewrite the total action in the ADM form 



S 



HP 



aSr + Sm — d x 



u'ijdtuj," + WjdtV^ + ndtX + udtX - {nh + N'^Ha + uj^'^g, 



IJr 



tIJ 



(13) 
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where the Hamiltonian, DifFeomorphism and the Gauss Constraints respectively read 

H = ^U'}j,U'^jLv''''Kb' - Dali'i + ^nlj {lia" DaX-D:ja"n) , (14) 

Ha = n'ijRj'' ^V^DbU'i + liDaX + DAIi, (15) 
Gtij = -DaIl'}j-^n^jVj]a + l(Xc7ijU-naijX). (16) 

Subsequently we will drop the subscript t from Gtu and denote it as Gij. From the definition of ajj and using the 
properties of gamma matrices (jA2p , Qjj can also be written as 

Gij = DaUlj + inf.l/j], + ^ (aA75%7J]A + euKLXl^N^'^l'^^x) « 0. (17) 

This is just a rewriting of the Gauss constraint in a form which will be useful in some later calculations. We 
have identified the basic variables in our theory. The gravitational degrees of freedom are incorporated in the pair 
(U°'jj,ujJ'''^ , which have 36 degrees of freedom and in the pair (llj,Va) , which have 24. The matter degrees of 
freedom are in the two pairs (ll, A) and (ll, A) each with 8 degrees of freedom. The total number of degrees of freedom 
without considering the constraints is therefore 76. 



B. Primary Constraints 

Let us now consider the constraints in the theory. At this stage we have the following constraints, (i) Since there 
is no momentum corresponding to w^^'^, we have to impose 6 primary constraints IIjj « 0. This leads to 6 secondary 
constraints Qjj w 0. (ii) Also there is no momentum corresponding to . We have to impose 4 constraints 11^ w 
which lead to 3 secondary constraints Ha ~ and 1 more secondary constraint H 0. (iii) From Eq. ([S]), we can get 
two other sets of primary constraints 

Cf := U''j-ae'''"'DbVci^O, (18) 
^ij - le^''eiJKLV,''V,^ « 0. (19) 

From (ITSl) we get 12 constraints, while ^T9\\ gives 18 because of the antisymmetry in IJ. (iv) From the definition of 
the momenta corresponding to the fermions (Eq. ([9])) we get 8 further constraints 

* := H-i^Nkj'' {l + i^l')x^O, 

* := TT + i^X (l + i^j^^ Nki''^ « 0. (20) 

These are the primary constraints of our theory. 

Including all of above constraints we can write the total Hamiltonian as 

Ht NH + N^-Ha + h'-'Gi.j + TaQ + A^'^$?j + m^- + (21) 

where the expressions of the constraints are given in equations (fT4|) . ([T5|) . (fTB]). (fTSl) . (fTO]) and ([20| respectively and 
the functions multiplying the constraints are the Lagrangian multipliers. At this point they are completely arbitrary. 
We now need to check whether the Hamiltonian system is consistent. To ensure the consistency of the Hamiltonian 
system, the constraints have to be preserved under evolution, i.e., the constraint algebra is closed. 

Note that the Gauss constraint Qu generates the 5*0(3, 1) transformations and is first class. However, as in [2l| 
the constraint which actually generates the spatial diffeomorphisms for the gravitational degrees of freedom is a 
combination given by 

Ha Ha + ioJ-'Gij + -eatcC^jWj (22) 
a 

This can be easily demonstrated as: 

S"-oj," {w/^ HaiN-)} = N-daUjJ-' + iVa'-'d.N- = CnW. 
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S"-nj := {nj, i/,(iV'')} = N'^dallj - W^daN' + U'jdaN'' = Cn^Wj. (23) 
For the matter degrees of freedom the constraint (1^^ acts as 

S"-X=\x,Ha{N'')}^N''daX , S"'Il= \ll,HaiNn} ^N^'dall + ndaN". (24) 

Clearly this combination Ha, acting on all the variables, generates Lie derivatives [l^l and can therefore be identified 
as the diffeomorphism constraint. Using the property of Lie derivatives (or by explicit calculation) it can be shown 
that H„ is first class. 



C. Second-Class Constraints and Dirac Brackets 

The terms in the constraint algebra which are not weakly zero are respectively 

{$?,(Ai^),i/(7V)} = i^M^.^NXj,VfrAyaXi' + euKLX^^), (25) 

{Q(7f),<fjA'(A^^)} - e^'^jlV^'iaXil' + ejjKLX^^), (26) 

{Cml), H{N)} = -^e'^'^^jiV^' (DA aijll - H ajjDA) , (27) 

{C!{ji)Mu)} = -^u A (l + z|75) ^J^/ (28) 

{cr(7a'), = r' (i + *f 7') x^i n?j, (29) 

{^(u),'^{u)} = 2iuuy/q-/^Ni. (30) 



As mentioned above tj/j and iJ^ are first-class constraints. We would also like the Hamiltonian constraint H to be 
first class. The remaining constraints are second class. Second-class constraints are problematic because the flows 
generated by them do not lie on the constraint surface. There are two ways of dealing with second-class constraints: 
First, after obtaining a consistent Hamiltonian system we can solve them and eliminate spurious degrees of freedom. 
Second, we can use Dirac Brackets instead of Poisson Brackets to define evolution and set the second-class constraints 
identically to be zero |23|. 

In this paper we shall construct Dirac Brackets only for the constraints 5* and while the remaining second-class 
constraints will be solved after obtaining a consistent Hamiltonian system. If we strictly followed the Dirac procedure, 
we would have to construct Dirac Brackets using all the second-class constraints. The algebra then became very 
cumbersome and complicated. On the other hand, Dirac Brackets for only ^I^ and ^' are easy to construct and, as 
we shall see later, greatly simplify the calculation. Recall that, for a set of second-class constraints {Xm}, the Dirac 
Bracket between two phase space functions / and g is given by 

{f,9r ■■= {f,9}~{f,X^}M-lAXn,9} (31) 

where the matrix Mm„ given by Mmn = {Xm,Xn} has a non-zero determinant. For the set of constraints 'I' and 
the 2x2 matrix M can be easily constructed from Eq (|30p . Subsequently, we shall use Dirac Brackets instead of 
Poisson Brackets and set the constraints ^ and ^ strongly equal to zero. Note that from definition ([31]) it is obvious 
that the Dirac Bracket is as same as Poisson Bracket if either f oi g commutes with 4' and ^ . 
Since the two constraints 4* and ^ have been eliminated, the total Hamiltonian (|2T]) becomes 

Ht NH + N'^Ha + h'-'Gij + jiCf + A^^^Jj. (32) 

Note that the evolutions generated by Poisson Brackets and Dirac Brackets differ only because of the 7aC° in Ht, 
since all other terms in Ht commute with \1/ and VE*. Also note that the gauge transformation and diffeomorphism 
generating characters of Qjj and Ha respectively do not change when we use Dirac Brackets, since they commute 
with * and W. 
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D. Secondary Constraints 



For a consistent Haniiltonian system, the constraints should be preserved under evolution, i.e., for all the constraints 
Cm, we require '■= {Cm, Ht}* ~ 0. Our analysis will be along the lines of [2l|. But owing to presence of fermions, 
it will turn out that torsion is not zero. As a consequence, the calculations are much more complicated, and the 
conditions of the closure of constraint algebra will now relate the torsion and the spin degrees of freedom. 

Let us first consider the constraint $77. Note that, since commutes with and in this case, its Dirac 
brackets are as same as its Poisson Brackets. Hence we have 

= {<^lj{ai'),H{N)] + {>fl,{a'J),C]U)] (33) 
where al^ is an arbitrary smearing function. Using Eqs ([SS]) and ([26|) . and after some calculation, Eq.(|33l) implies 

Hence we have 

HV'' - llV^ - livl + livi) - eafcc^ (iV^n''-^ - ^^H'^^ - ^A75 [V^'l' - V^' l'\ a) « 0. (34) 
Multiplying with and using the properties (jBSp we get 

27c + yhiyl - Vhi + eabc-N^n-'vH - eabc^Xl5 {V^^Vj^^ - V'-'VjW) A « 0. (35) 

a Z 

By multiplying this equation with Nj, Vf and Vj respectively and using the relations (|B2p and (|B3p . we obtain the 
following relations 

jiNi = ^eabcU'}Vj\ (36) 
Za 

liVf = 0, (37) 
liVdi = edbc^Xl5ViV\ (38) 

where we have used Eg. ([57)1 to obtain Eq. ([55)) . Finally from the equations and ((55)) we get a solution for the 
Lagrangian multiplier 7^ as 

7c' = ^abc^Xl^V^Vh'X - ^eabcN'll'}V,l (39) 

Note that, all these equations differ from the corresponding equations in [2]| only by the fermion-dependent terms. 

So, we have obtained 12 components of 7^ from the 18 equations in (|34)). Consequently there are 6 constraints 
remaining. By inserting the solutions ((5^ back into Eg. dMl) and after some calculation, we get the following secondary 
constraint: 



^ab _ ^ayb ^ ^bya _ ^^X^^V^V} N Kl'' \ ~ 0. (40) 

Since x"*" is symmetric in (a -f-)- b), it contains just the 6 required constraints. Before proceeding further, let us 
count the degrees of freedom. As mentioned above, the total number of degrees of freedom without considering the 
constraints is 76. Now let us count the degrees of freedom removed by the constraints. Clearly Qjj and Ha are first 
class. We would like H to also be first class which would give us (6 + 3 + 1) = 10 first-class constraints removing 20 
degrees of freedom. The constraints Cf and j are primary second class removing 12 + 18 = 30 degrees of freedom. 
By setting 5* and ^' strongly equal to zero we have removed 8 more degrees of freedom. Finally the secondary 
constraint x""^ turns out to be second class and thus removes 6 degrees of freedom. Thus the number of independent 
degrees of freedom in our system is 12. 



8 



E. Relations and Simplifications 



We have found all the constraints in our theory. We now need to show that the Hamiltonian system is consistent, 
i.e., the constraints are preserved under evolution by H^, where the evolution is defined by using the Dirac Brackets 
which we have constructed. This will put restrictions on the Lagrange multipliers. We shall solve the Lagrange 
multipliers of the second-class constraints and keep those of Qij, Ha and H free. As seen above, $jj ~ fixed the 
Lagrange multipliers 7^ of the constraint Cj to the form given by Eq. p9p . This can however be further simplified. 
For this and for subsequent calculations, we now derive some useful identities using the constraint equations. All 
these identities hold weakly, i.e., they are true only when the constraints are used. 

From the constraints (|18p and ([T9l) we can easily obtain the relation: 

Dij := DJilj - -euKL^-'^'^V^ « . (41) 
a 

Using this and the Gaussian constraint pT)) we get, after some algebra, 

nfz^j]a + ^A75iV[/7J]A « 0. (42) 
Multiplying this equation with Nj and then with Vj and using the properties (|B2|) and (|B3[) . we get 



^'jN' « '^Xl.ljV^'X. (43) 
By multiplying relation (|^^ with and then with Vj and again using the properties (IB2[) and (|B3[) . we get 

njt//-njv7 w 0. (44) 

Using Eq. (^H) in the constraint (|40l) we get the relation 

n?V/ « '^\l5Vtv!NKl''\. (45) 

These identities can be used to greatly simplify the subsequent calculations. 

First, note that because of the identity (j44|) . the second term on the RHS of Eq. (p9|) drops out and the Lagrangian 
multiplier of Cj in Ht becomes 

li - eabc^Xl.VtV'j^'X. (46) 

This leads to further simplification of our problem. As mentioned above, the difference in evolutions generated by 
Poisson and Dirac Brackets comes only from the 7^Cj in Ht, where 7^ is an arbitrary Lagrangian multiplier. The 
condition that $f j be preserved under evolution has fixed 7^ to the specific form given by Eq. p6)) . Now recall from 
Eq. ()28p . for an arbitrary smearing function 77^ we have 

When f]^ = 7a, which is of the form given in Eq. using (|B5[) and (|B2[) we get 

li ^ij = (ea,c^Xj,j''X^ VfVj^V(jNj] = 0. 

Using this we can see that {C" (7^), \E'(m)} ~ 0. Similarly, we also have {Cj (7^), 4'(m)} « 0. Therefore, once 
the Lagrange multiplier 7^ is fixed to the value required for a consistent Hamiltonian system, we can calculate the 
evolution using Poisson Brackets, i.e., for any function / on the phase space we have 

/ = {LHtV « {LHt}. (47) 

This key simplification further justifies our choice of Dirac Brackets. 

Second, consider the identity (j45|) again. Multiplying it by Vbi and using the properties (|B3P and (jBSP , we get 

n? « ^A75 (V[?iV,,]) 7^ A = |A75n?^7-'A. (48) 
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This equation relates the torsion degrees of freedom encoded in 11" with the spin degrees of freedom A and A. Note 
that we have used only constraint equations and not equations of motion in deriving Eq. ()48|) . This is a weak relation 
which comes naturally as a condition for the closure of the constraint algebra, since it has been derived by using the 
secondary constraint x""^- When there is no matter, this equation would indicate that torsion is zero [2l|. Note that 
relation p8)) is same as the one obtained in [Tsj . 



F. Closure of Constraint Algebra 

At this point let us see what are the remaining requirements for obtaining a consistent Hamiltonian system. Recall 
that we had started with a constraint algebra with non-zero terms given in Eqs. (|25II30|) . However, instead of Poisson 
Brackets we used Dirac Brackets constructed by using ^I^ and 5'. Since H and $fj commutes with and 5*, the 
constraint algebra terms (|25ll27p remain the same. Also, now that = and = are strong equalities, we do not 
need to worry about terms (|28ll30p . Ensuring that the RHS of Eqs. (|25ll26p are weakly zero gave us a new constraint 
X"''' and fixed the Lagrangian multiplier 7^. The new secondary constraint turns out to be second class, and we 
have additional non-zero terms in the constraint algebra as 

{^hi^c').x''\^ab)} = 2aabX'/e-'''eijKLVl^V,'^, (49) 
{x"''Kfc),*(^)} = 2iuaabVj''V''N,n'X, (50) 
{x^''{aab),^{u} = -2iuaabVfV''NjW, (51) 

{x''H<yab),C^j{li)} - ^e^H'Vb^ (nf/iV,/] - aUljXj.NKl^x) (52) 

+^iiNj (2n^y;n^, - ^n'},v^u'},Xj,j^x) - ^e-''7v,nj,i?.7i, 

D, (^^nf,|^|nj]^) - ^A^[,n}] (53) 



{x^\aab),H{N)} = ^ n^jNj^^aU]jXj,NMl^'X 



2^9 



2q 

{m'i.NiVj^U'kj - an}in^^iF^iV,A757''A) {NN') 



2^/9 



Again, since we are using Dirac Brackets we do not consider Eqs. ([50]) and ([5T|l . Also, as mentioned above, the 
Poisson Brackets are weakly equal to the Dirac Brackets. The RHS of remaining equations have to be made weakly 
zero. The subsequent calculation is long and complicated. Here we only describe the procedure and more details are 
provided in the Appendix ([C|) . 

Since is a secondary second-class constraint, we do not add it to the total Hamiltonian Ht which is still given 
by Eq. Now in Eq. ([32]), iV, 7V°, A'^"' and A^"' are Lagrange multipliers which can be arbitrary functions ,but 

7^ is fixed to the specific form (pS)) to ensure sa 0. Now we need to ensure 

Cnvi) = {cnvi),HT} « {cnvi),{^ljiXi') + HiN))}^0, (54) 

r'iTab) = {x"''(^afc),ffT} « {X"''(^afc), {^IjiXi' ) + Cf (ji) + H {N))} « 0, (55) 

where rj^ and aab are arbitrary smearing functions. We solve the 18 independent equations ()54p and ()55p to fix the 
18 independent components of the Lagrangian multiplier A^"^ (see Appendix jCj). Note that the constraint H is still 
second class. However it can be shown that we can construct a first-class Hamiltonian constraint by 

H = H+:^Cf + ^^'}j. (56) 

This will go back to H once we solve the second-class constraints in the next section. 

It can be explicitly checked that the constraint algebra is now closed. We have obtained a consistent Hamiltonian 
system. The Lagrangian multipliers of the second-class constraints Cf and 4>Jj have been fixed by the conditions of 
closure of the constraint algebra, while those of the first-class constraints Sij, Ha and H remain free. 
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IV. SOLVING THE SECOND-CLASS CONSTRAINTS 



Having obtained a consistent Hamiltonian system we now proceed to solve all the second-class constraints , <& " j 
and x"'''- We will do this after performing a partial gauge fixing. Since we have already proved the consistency of the 
Hamiltonian system, we can be sure that making a gauge choice now will not lead to any inconsistency. Our goal 
is to reduce the internal SO{3, 1) gauge symmetry to SU{2). So we break the S'0(3, 1) by fixing the internal vector 
N-^ = (1,0,0,0), i.e., we fix a specific timelike direction in the internal space. This is a standard gauge choice and is 
known as time gauge. From Eqs. (jBl|) . (jB2|) and (|B3| it is easy to see that 



TV-f = (1,0,0,0) ^ 1/0 = = V^. (57) 
For consistency, this gauge fixing condition has to be preserved, i.e.. 

Hence in time gauge we get 

A°V„ « daN. (58) 

The Lagrangian multiplier of Qoi gets fixed. This is expected because, by fixing , we have broken the 50(3, 1) 
gauge invariance. The preservation of this gauge fixing condition implies that the boost part of the Gaussian constraint 
does not generate gauge transformations. 

We first solve the constraint p^ . which can be written as 



Thus in time gauge we have 



n^,— ; = \e'^'^e,,uV^V^ -.^Ef. (59) 



So, after solving this constraint only the Hg^ part of H°j remains a dynamical variable. Consequently, only the uj^^ 
part of the S'0(3, 1) connection remains dynamical. 

For convenience we define if* := 2w^0' which will be conjugate to Ef . The uj^^ is the remaining part of the 
connection which will get solved in terms of other variables while solving the remaining constraints. Since our gauge 
group is now reduced to SU{2), we will expand 5*0(3, 1) connection components uj^^ in the adjoint basis of SU{2) 
as uj^^ :— —e^^'^Tak- The quantity Tak is the SU{2) connection. Note that we had started with a S'0(3, 1) spin 
connection uj^^^ which was not torsion-free. Therefore the variables and FJ^ that we define above will contain 
information about torsion implicitly. Also, from (Ilj3|) it is clear that, in time gauge, is the inverse of V^. Using 
the properties of inverses and determinants of matrices, it is easy to see from ([5^ that Ef = ^/qVf- is the densitized 
triad. We can also determine its inverse E^ = ^V„\ ^ 



Next we consider the constraint (|40|) . Actually, we have already found the solution of this constraint by Eq. 
It can be verified that the relation (j48| identically solves the constraint. Using the above solution ([59| . in time gauge 
the solution P5|) simplifies to 

ng-fi^rWA ; H," = -|£;rA757°A. (60) 

The torsion degrees of freedom are solved in terms of the densitized triad E^ and the fermionic fields A and A. 
Using above results we can now solve the remaining second-class constraint ([18]) as 

CS = f^(wA--^e.,feX^i?''-) -0, (61) 

Of = ^ ^A757°A + -i= [e^^'E'^ElEtd.El + \e^^^ E^ ElE%Ei^ + {tIeIE^ - TIe\EI^ = 0. 

(62) 



^ In this work, we are not interested in the behaviour under parity transformations. Therefore we omit the sgn(det_E?) terms from our 
expressions. 
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Equation (|62]) can be used to solve the spin connection FJ^ in terms of the other variables. After some algebra we 
obtain 

= \^"^EiElE'id,Ei + \e^'''E',E]Eld,E[ + ^e^^^Eld^El - ]^e^=^Eld,Ei - ^ElX^.^'^X (63) 
K-^El\^,f\ (64) 

where we have denoted by F^, the first four terms of which do not depend on the fermions. It turns out that F^ 
is exactly the SU{2) spin connection which we would have obtained had there been no fermionic matter 5]. So, when 
there is no matter we go back to the standard formulation. Also note that the spin connection F^ is independent of 
the arbitrary coupling parameter a. 

Let us now consider the Gaussian constraint pT|) . It is obvious that because of the identity the Gaussian 

constraint is given by 

Gij = Dall'ij + ^ (e7jKLA75AfI^7'^'A) « 0. (65) 
In time gauge and using the solutions of ^ we can rewrite this as 

Go; - daE'^ + Lj:^E;;^daEf + e,jkriE'^''-0, (66) 
Gjk = ^K.^E',^ - yle.kiXWX « 0. (67) 

On rewriting oj^^^ in terms of FJ^ and inserting the solution (1551) . the Goi term becomes identically zero, i.e., DaEf — 0. 
The boost part of the Gaussian constraint, which has no physical meaning after choosing time gauge, is therefore 
eliminated from our theory. Also comparing Eq. (|67p with Eq. (|5T|) it is easy to see that if ejkiGjk ~ then Cq is 
automatically zero (assuming ^ 0). 

So far, we have reduced our original phase space by consistently imposing time gauge and then solving all the 
second-class constraints. As a result, some basic variables in the original phase space have been eliminated in terms 
of the others. To obtain the basic variables in this phase space we need to study the pq term in Eq. (|13l) . 

nijdtujj' + Wjdt V/ +TldtX + UdtX. (68) 

Using the solution ([591 the first term becomes EfdtK^. Using the solution (|60l) and the fact the preservation of time 
gauge requires that dtV^ = 0, the second term becomes —^X^j^XE^-dtV^. For the last two terms in (|68p . recall that 
A := A^7°. Also 11 and 11 can be read off from the identities ([20]) . Then in time gauge, using the properties of the 7 
matrices given in appendix (fA| we get 

na^A + ndtX = -i^x^ {dtX) + iVq{dtX^) a + ^x^^x 9* (Vq) ~ Ot (|V9A^75A 

Now since 

dtVq = V^V,'^dtV: = EtdtV:, 

we have 

n^tA + na* A = -i^/^x^ (dtX) + (dtX^) x + ^x^-f^XE^dtV:, (69) 

where we have neglected the total time derivative term. Putting everything together, (|68p becomes 

E^dtK - *V9A^ {dtX) + {dtX^) X 

= EtdtK + (5tC^) K - [dtO 

^ EtdtK + {dtC)li^,+IiddtO, (70) 

where, following [i^IIIIj have defined half-densities of the fermionic variables: C, := .^A and '■— -^A^. In 
the third line we have identified 11^^ = ^^C^ ^ind Il^t — iC,- The basic variables of our system after solving all the 
second-class constraints are clearly shown in the expression (j70p . By using the half-densities, we have removed the 
,Jq factors from the fermionic terms. Consequently, the fermionic sector and the gravitational sector now commute 
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with each other triviahy. Therefore the Dirac brackets can only differ from the Poisson Brackets for the fermionic 
variables. The constraints ^ and which we had used to construct the Dirac Brackets now become 

V' := + iC^ « ; ?/i B^t - < « 0. 

They satisfy 

{^,V^} = 2i. (71) 

So the determinant of the required matrix Mmn (see (|3ip ) is not zero. Using this fact we can calculate the Dirac 
brackets for the fermion sector. The remaining 7 constraints, Qjk, Ha and H, are first class. 



V. SU{2) GAUGE THEORY 

We have obtained a consistent Bamiltonian system which is invariant under local SU (2) rotations. Bowever this 
is not a SU{2) gauge theory yet. The basic variables in the gravitational sector are the densitized triad Ef and its 
conjugate K^. The connection F^, given by Eq. (|63p is a function of Ef, C and C^. In this section, we shall make a 
canonical transformation to find a new variable conjugate to E°' and which also turns out to be a connection. 

To do so, let us look at the Gaussian constraints (|66| and (|67p again. Recall that the first equation is identically 
zero while the second equation is a first-class constraint. We construct a new combination 

GQ^ + iSe.jkG'" 
= daEf + e„fe in + pKi) E'^' - /3V^A757*A « , 

where (3 is any non-zero real number. Using Eq. (j64p we get the standard SU{2) Gaussian constraint with matter as 

= - /3CH%57»C « 0, (72) 

where :=T\-\- . Tensorially, the new connection which we have defined is in the same form as the standard 
Ashtekar-Barbero connection without torsion. The basic variable B° which encoded the torsion has been solved in 
terms of the fermionic degrees of freedom when we solved the second-class constraint x"'^ in Eq. (pO)) . We had started 
with a 5*0(3, 1) connection lu^'^ which is not torsion free. That fact is refiected in our expression of the SU{2) spin 
connection in Eq. (j64p . But in the new connection which we define above, we remove exactly that additional 
piece. Bowever, since we had defined := 2a;„°*, the variable implicitly contains information about the torsion. 
When there is no matter, torsion goes to zero and the St term in our action ([T]), and therefore, the terms originating 
from it in the Bamiltonian analysis vanish j2ll |. Then we go back to the standard formalism with a torsion- free 
50(3, 1) spin connection. 

With our connection A^, we can also easily prove the canonical commutation relations following [5] (with k = 1) as 

{Enx),Al{y)} ^ 6tSiS{x,y), (73) 

{Enx),E^{y)} = = {A:,{x),Al{y)}. 

We have obtained a SU{2) gauge theory formulation of our system. A peculiar result of our analysis is that, the 
coupling constant a of St in the action ([1]) is absent from the definition of Al^. It depends on an arbitrary number /3 
which may not have any relation with the parameter a with which we started. Bowever, in the analysis performed 
with Bolst term or the Nieh-Yan term , the connection depends on the arbitrary coupling parameter 7 present in the 
action. ^ 

The remaining constraints can also be rewritten in terms of the new basic variables. Using — ^{A^ — F^) and 
the Gaussian constraint ([7^. the diffeomorphism constraint ([T5|) can be written as 

Ha = E^d[aKl]~KdtE'^+^(^TT^daC-ida7Tc)C+{daC^)7rcf-C^dan(f 

« ^E'^F:, - ALC^7%57.C + I (^HdaC - (dancK + (5aC^)^ct - CdaTT^f^ , (74) 
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where := + ejj.A{A^ is the curvature of A^. This is exactly the standard diffeomorphism constraint. Note 

that, although the arbitrary parameter /3 appears, the a dependent terms have again dropped off from the final 
expression. The Hamiltonian constraint ([Ti| is more complicated. After some calculation, we get 



1 



+ — ( ^cTsC - C^757rct ) e^.kEfEpaE, 



(75) 



where Ti 



\ai and Ui are Pauli matrices. Again this expression is independent of a. This expression goes over to 



the standard expression when the fcrmions are set to zero. Thus we complete our task of obtaining a SU{2) gauge 
theory. 



VI. CONCLUSION 



Let us briefly summarize what we have achieved in this paper. We started with an action containing a torsion- 
squared term and fermionic matter apart from the standard Hilbert-Palatini term. This term is just the difference 
between the total derivative Nieh-Yan term and the Hoist term. Since an SU{2) gauge theory formulation can be 
derived from actions containing either [l^, [l^ , it seemed possible that such a formulation can also be obtained from 
our action containing only the term. We also need to add fermionic matter because the vacuum case is torsion free 
[2]| and we are left with only the Hilbert-Palatini part which docs not admit a gauge theory formulation. We take 
non-minimally coupled fermionic matter so that the classical equations of motion for the fermions do not depend on 
the coupling constant a multiplying the torsion term. 

We do a 3 -|- 1 decomposition of our action, do a constraint analysis and finally obtain a consistent Hamiltonian 
system with second-class constraints. The second-class constraints related to the definition of fermion momentum 
are treated in the standard way via Dirac Brackets, while the remaining second-class constraints are solved after 
breaking the SO{3, 1) invariance by fixing time gauge. As far as we know, such Hamiltonian analysis on an action 
with non-zero torsion term with explicit expressions of all the second-class constraints is new in literature. Similar 
analysis with the Hoist term (with non-minimally coupled fermions) [17rrl9j and the Nieh-Yan term (with minimally 
coupled fermions) 15] has already been attempted before, although the picture is not yet complete. Apart from the 
crucial fact that the gravitational part of our action being different from those studied in literature so far, there are 
several other differences in our approach. Since we are motivated by PGT where the initial action is invariant under 
Poincare transformations, our starting variables are different from those used in [13) [HI- Unlike the treatment in 
[To! we do not break up our variables into the torsion dependent and independent pieces. Moreover, since we do not 



have the Hoist term, the techniques developed in 12211 for dealing with second-class constraints and used in [15[ are 
not available to us. Also, unlike the treatment in we fix the time gauge only after we have found all the 

second-class constraints and obtained a consistent Hamiltonian system. Also, in [l9[, there are two different coupling 
parameters, one for the Hoist term 7 which is different in general from the parameter a in the non-minimal fermion 
action. In our case the parameter in the term gets related to the one in the fermion action due to the closure of 
constraint algebra. 

On solving the second-class constraint x'^'', torsion gets related to the fermionic degrees of freedom via Eq. pS)) 
which is as same as the one obtained in ^5|. Further, solution of the second-class constraint Cf gives the SU{2) 
spin connection in terms of the densitized triad. This differs from the spin connection in GR ^ only by a term 
which depends on the fermions. In the final step we obtain the connection dynamics by defining a new connection 
Al^ which is algebraically in the same form as the Ashtekar-Barbero connection without torsion. However, unlike the 
torsion-free case, the part comes from the ujj^'^ part of the 5*0(3, 1) connection which is not torsion free. As a result 
it is not obvious the K^^ can be directly related to the extrinsic curvature Kab on shell. While the diffeomorphism 
constraint (|74p is standard, the additional terms in our Hamiltonian constraint (|75p are somehow different from the 
ones obtained in literature. Although these constraints can be loop quantized using existing techniques, it may be 
possible to rewrite them in a form more convenient for loop quantization. We leave this issue for future research. This 
present work at least opens the door to extending loop quantization techniques from standard GR to more general 
PGT of gravity 

One more peculiar and novel feature of our analysis is that the coupling constant a in our starting action totally 
disappears from the final Hamiltonian system. Rather the new connection depends on an arbitrary new parameter 
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/3 which does not have any relation with a. This is unhke all the analysis done with Hoist and Nieh-Yan terms 
where the Immirzi parameter 7 in the definition of the Ashtekar-Barbero connection is just the coupling parameter 
in the action (up to numerical factors). The nature and the role of this new parameter is not clear so far and deserve 
investigating to other actions with torsion terms. 
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Appendix A: Gamma Matrix 

In this section we collect some of the standard properties of Dirac matrices which we have used in previous sections. 
The 7 matrices, in any dimension, satisfy the Clifford algebra 

{ill 7j} = 7/7J + 7J7/ = 277/j (Al) 

where ?7/j is the flat Minkowski metric. We shall restrict ourselves to 4 dimensions and choose the signature ( — (- ++) 
which is different from the signature usually used in QFT. In this signature the above relation can be decomposed as 

70' = -I4 ; 7.' = I4. 

This implies that 70 is anti-Hermitian while 7^ is Hermitian. Note that all the 7 matrices are unitary. We also define 
the commutator cr/j ^[7/, 7j] and another standard combination 75 :— i7o7i7273- It is easy to check that (75)^^ = I 
and (75)^ = 75. In the Weyl representation, commonly used for massless fermions, the Dirac matrices can be explicitly 
written as 

il, \ / -ioi \ / -I, \ 

j ' j ; ^^=(, l2j- 

In this paper we have used the following standard identities 

{75,7/} = = [75, CT/,/]; 

\lK,oij\^nK{ilj\ ; {7x,cr/j} = ie^^'^^757L. (A2) 



Appendix B: 3 + 1 Decomposition 

In this section we give the parametrization of the tetrad and the co-tetrad fields which we use in this paper. They 
read 

eti = NNj + N'^Vai ; e" = - — , 

eai^Vai ; 6"^ = ^^'^^ + ^^, (Bl) 
with N^Vai = ; N'Ni = -1. (B2) 



What we have done is that we have reparametrized the 16 degrees of freedom of e^/ into 20 fields given by (jBl[) 
subject to the 4 constraints (|B2[) . From these definitions, the following identities also hold: 

V'Vu^S^ ; V^"'iV/ = ; Na-.^VaiV.'N", 

y'Va = v''' + N'N''. (B3) 
In terms of these fields the metric takes the standard form 

_ / -N"^ + N^Na Na 



15 



It is easy to see that 

g det(5^,) = - dci{VaiVl) , 

e := |det(e^/)| = N ^ dei{VaiV^') = N^det{qab) = N^. 
Using the definitions given above we can also prove the fohowing two identities which have been used in our analysis, 



tt[i KL , AT[o.Tf] 



/J' 



(B4) 
(B5) 



Appendix C: Determination of A^'' 

In this section we show how to obtain A^'-' from Eqs. ((54)) and fSS]) . First let us consider Cf. Using Eqs. ([26)) and 
27)) . we get 



aN 



For convenience, we define 



aN 



Xy := {aX'J + euKLX^"-) - — {DaX aijU - H aijDaX) . 



(CI) 



(C2) 



This is the quantity we shall solve for. The equation ()C2I) can easily be inverted to express A^^^ in terms of X^'' as 
X" 



^ {DaX aijU - n aijDaX) - ^cijkl {DaX a'^^^Il - U a'^'^DaX) + X'J - -euKLX^"" 



a2 + 4 



Thanks to this definition, it is easy to see from (|C1[) that 
Now let us consider x""^ . We have 

r'K^) - {x''\c^ab),^h{>^a)] + {x^\<ya,).{Cnii)+H{N))] 



= —'^CTabXl'^ e°''^'^eijKLVj^Vj" + aab'S"'" , ^ era, 



b-irL 



(C3) 



(C4) 



(C5) 



where we have defined (Tab^""^ '■= {x"''(fafc), (C/(7i) + ^(-^))}- The explicit form of E"'' is very complicated and 
can be calculated using Eqs. (|52p and (|53)) . However we do not need the explicit form of E"''. We are interested in 
solving for A^'' which only comes from the first part of Eq. ()C5I) . After some more algebra we obtain the equation in 
terms of Xl'' as 



a 



+4 



2Vfvfn'k^x^^ - (vrifn^i + v^v^n^A x'. 



^cab f^^ 



-e-" y—eijKL {DaX a'^^n - H a'^'^DaX) + — {DaX a'-'li - H a'^DaX) ] VfV^ 



4N 



-e"^"" (^^iJKL {DaX a^"n - n a^^Z?aA) + ^ {DaX a^'^H - H a'^DaX) ) F/T// 
« 0. 



(C6) 
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Using (jC4l) and (|C6|) . after a long calculation we get 
1 



X 



ij 



4V9 



4a^ 



(C7) 



where, for brevity of notation, we have defined 



Ai' ^eijKL (DA a^^^H - H a^^D.A) + ^ (D:A a^-^H - H a^'^i^.A) 



Putting dCTl) into ([Cll, we get A^"'. 
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